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with positive kinetic energy (coupled quintessence) and with negative kinetic en- 
ergy (phantom). In the two cases, we investigate the existence of attractor solu- 
tions which correspond to an equation of state parameter uj — —1 and a cosmic 
density parameter fl^ = 1. We find that the coupled term between matter and 
dilaton can't affect the existence of attractor solutions. In the Mexican hat po- 
tential, the attractor behaviors, the evolution of state parameter uj and cosmic 
0^ ' density parameter fi, are shown mathematically. Finally, we show the effect of 

coupling t 
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1. Introduction 



■ , Observational data from SNe Ia[l], WMAP[2], SDSS[3], show that we are living in a spatially flat universe 

i^' which consists of about two thirds dark energy with negative pressure, one third dust matter including cold 

(~| , dark matters plus baryons, and negligible radiation, and that our universe is undergoing an accelerated 

expansion. In order to explain current acceleration of the universe, we require an exotic energy dubbed 
" dark energy" with equation of state satisfying a; < — | . 
^^ ■ Many models have been proposed so far to fit these observations. Among these models, the most typical 

^ , ones are cosmological constant and a time varying scalar field with positive or negative kinetic energy evolving 

in a specific potential, referred to as "quintessence" [4-14] or "phantom" [15-19]. The essential characteristics 
of these dark energy models are contained in the parameter of its equation of state, p = tup, where p and p 
denote the pressure and energy density of dark energy, respectively, and w is a state parameter. Quintessence 
model has been widely studied, and its state parameter u which is time-dependent, is greater than —1. At 
the same time, observations also suggest another dark energy model — phantom with an equation of state 
a; < — 1. The most typical characteristic of the phantom model is that the kinetic of the scalar field is 
negative. Because of this, it leads to some unusual features of phantom model. A striking consequence 
of phantom model is that the Universe will undergo a catastrophic "Big Rip" in a finite time. However, 
Hao and Li [20] have presented that a de Sitter attractor will prevent the phantom energy from increasing 
up to infinite in a finite cosmic time, therefor the presence of phantom energy dose not lead to a cosmic 
doomsday in a theory with de Sitter attractor at late time. On the other hand, it has been shown that big 
rip singularity can be generically avoided if one considers quantum gravitational corrections to dynamics [21]. 
Some authors have considered the coupled dark energy [22] . In the work of L. Amendola[23] , he investigated 
the cosmological consequences of the coupled quintessence model, assuming an exponential potential. In 
our previous paper [24], we have considered a dilatonic dark energy model, based on Weyl-scaled induced 
gravitational theory. In that paper, we find that when the dilaton field is not gravitational clustered at 
small scales, the effect of dilaton can not change the evolutionary law of baryon density perturbation, and 
the density perturbation can grow from z ~ 10^ to z ^ 5, which guarantees the structure formation. When 
dilaton energy is very small compared the matter energy, potential energy of dilaton field can be neglected. 
In this case, the solution of cosmological scale a has been found[25]. In what follows, we shall consider 
a more general dilatonic dark energy model that the contributions from radiation and matter can't be 
neglected. In this model, we will study coupled quintessence — dilaton field with positive kinetic energy, and 



phantom — dilaton field with negative kinetic energy and find out the sufiicient conditions of existence of 
attractor solutions, which limit the choice of potential of scalar field. When we take the potential as the 
form ^(<T^ — e^)^ + Wq named Mexican hat potential, we investigate existence and dependence on initial 
conditions for a late time attractor and obtain the sufficient conditions for the existence of an attractor 
solutions. We find that these attractor solutions correspond to an equation of state w = — 1 and a cosmic 
density parameter 17^ = 1, which are important features for a dark energy model that can meet the current 
observations. These results are shown mathematically. Finally, we show the effect of coupling term on the 
evolution of X and Y with respect to N numerically. 

This paper is organized as follows: In section 2, wc introduce Weyl-scaled induced gravitational theory 
and regard dilaton as quintessence. Based on this, we investigate the sufficient condition for existence of 
attractor solution. In section 3, we consider phantom model and give the sufficient condition for existence 
of attractor solution. Section 4 is a special example — Mexican hat potential. Section 5 is summery. 

2. Coupled Quintessence 

To deduce the field equation of induced gravitational theory, let us consider the action of Jordan-Brans- 
Dicke theory firstly 

s^ J d^x^/^m - u-f^-M^ - A{cj)) + Lflu^dm (1) 

where the lagrangian density of cosmic fiuid Lfiuidi^) = ■^j'^'^df^ipd,,ip — l^(V')j 7 is the determinant of 7^,^ 
which is Jordan metric, uj is the dimensionless coupling parameter, R is the contracted Rp,i,. The metric sign 
convention is(-, +,+,+). The quantity A(0) is a nontrivial potential of (j) field. When A(0) 7^ the action of 
Eq.(l) describes the induced gravity. The energy density of cosmic fluid p = ^{-&)'^ + V{ip) and the pressure 

However it is often useful to write the action in terms of the conformally related Einstein metric. We 
introduce the dilaton field a and conformal transformation as follows 

<P = le"- (2) 

7m- - e-"-g^, (3) 

where a^ — „ "^ „ with w > 3500 [26] being an important parameter in Weyl-scaled induced gravitational 
theory, 5^^ is the Pauli metric. In this paper, we work in units(«;^ = 8irG = c = 1). From the solar system 
tests, the current constrain is a^ < 0.001 [27]. The new constrain on the parameter is a^ < 0.0001 [28], which 
seems to argue against the existence of long-range scalars. Perhaps such a pessimistic interpretation of the 
limit is premature [27]. 

The action(l) becomes Eq.(4) by performing the conformal transformation Eq.(2) and Eq.(3) 

S = J d^X^^[^R{g^,) - ^g^-d^ad,a - W{,j) + Lfiu^M)] (4) 

where Lfiuidii^) — ■^g'^^e^"'^df^^pd,yip — e^^°"^V{ip). The conventional Einstein gravity limit occurs as a — > 
for an arbitrary lu or uj ^ 00 with an arbitrary a. 

The nontrivial potential of the a field, W{(t) can be a metric scale form of A(0). Otherwise, one can start 
from Eq.(4), and define W{<7) as an arbitrary nontrivial potential, g^i, is the pauli metric. Damour and Cho 
et.al pointed out that the pauli metric can represent the massless spin-two graviton in induced gravitational 
theory. Cho also pointed out that in the compactification of Kaluza-Klein theory, the physical metric must 
be identified as the pauli metric because of the the wrong sign of the kinetic energy term of the scalar field 
in the Jordan frame [29]. The dilaton field appears in string theory naturally. 

By varying the action Eq.(4), one can obtain the field equations of Weyl-scaled induced gravitational 

theory. 
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+ e-"'^[a^V5.V' ' \9i.vdpWiA - g^.e-^^^FCV')} (5) 

Aa = ^^^ - ^e-2"-gA<-a^^9,7A - 2ae-2-T/(^) (6) 



dip 



AV = -ag^,d^i;d,<J + e'^'^ -^ (7) 



The energy-momentum tensor T^^ of cosmic fluid is 

Tt_,i, = {p + p)Ufj^U,y + pg^^ (8) 

where the density of energy 

P=^^'+e-""F(V') (9) 

the pressure 

P=\^^-e-^''V{i:) (10) 

p and p are related to their directly measurable counterparts by p = e^"°'p, p — e^^^p. 

In FRW universe, the field equations of Weyl-scaled induced gravitational theory can be expressed as 
follows: 

H'^l[la' + W{a)+e-'^^p] (11) 

a + 3HcT+^ = l-ae~°"'{p-3p) (12) 

d(T 2 

p + 3H{p + p) = ^a&{p + 3p) (13) 

where H is Hubble parameter. From Eq.(12), we note that there a coupling between matter and dilaton. 
Its effect will be discussed at the end of this section. For radiation pr = 3pr, we get pr oc ^^ from Eq.(13). 

For matter pm ~ 0, we get pm oc ^-3— from Eq.(13). Taking these results into Eqs.(ll) and (12), wc obtain 

H^ = Hfi^^l^^^ + n^^.e-^-C-^f + n^4"-r] (14) 

Pc.i a a 

a + 3H&+-— = -aT]ie-'!°"'a-^ (15) 

dcr 2 

where Hf — ^^, pd is the critical energy density of the universe at initial time ti. Hi, ^m,i, ^r,i denote the 

Hubble parameter, matter energy density parameter, radiation energy density parameter at initial time U 

3 
respectively, and rji = ^7''°' with pm i being matter energy density at initial time ti. We define our starting 

point as the equipartition epoch, at which Qm,i — ^r,i — 0.5 and consider the initial scale factor a^ — 1 for 
convenience. So we have 

^2 ^ ^2[M + ]^ ^ r!„^,e-i-a-3 + r!.,.a-4] (16) 

Pc.i 

The effective density p^r and effective pressure p^r can be expressed as follows 

Pa = \o''^W{a) (17) 

Pa = \o'' - W{o) (18) 

So, the equation of state of dilaton field is 

In order to gain more insights into the dynamical system, we introduce the new dimcnsionlcss variables 

X = — , r=-^, N^lna (20) 



\a^ - W{a) 



Eq.(19) becomes 
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WiX) 



^a = -^^ ^ (21) 



The field Eqs.(14)(15) could be rewritten as follows 

dX _ croF 



(22) 



dN^ ^^ ^3^ + AH ^^^> 

where the prime denotes the derivative with respect to X and H could be expressed as 

H = H^^ ¥tY' + W{X) ^ ^^^^^-i..„x^-3N ^ j^^^^„-4^]^ (24) 

Pc.i 

According to reference [30], the dynamic system decided by Eqs.(22)(23) belongs to the class of so called 
autonomous systems. The coordinates of critical point is decided by the following condition 



In dilaton-dominant epoch, the matter energy density can be neglected comparing to the dilaton energy 

4:H~ 



density, that is ^^^^^^ — jjj > 0. So autonomous systems (22) (23) become 





dX 


ctqY 






dN 


H 




dY 

dN 


^-3Y 


W 

0", 


'{X) 



(26) 

(27) 

The critical point of the above autonomous system is (Xc, 0), where Xc is defined by W'{Xc) = 0. Linearize 
the above equations about the critical point, we have 

d_X ^ V3a,Y 

dN ^/WiX:) 

^ = _3r - ^^"(^'^)^ (29) 

dN al^Wix:) 

The eigenvalues of the system are 



2^1,2 = ^ (30) 

where ^ = 3 and ( ~ WTxia^- ^'^^ ^ positive potential, if W"{Xc) > 0, i.e. xi < and X2 < 0, the critical 
point is a stable point which corresponds to a late time attractor solution. So, W{Xc) > 0, W'{Xc) — 
and W"{Xc) > is the sufficient condition for existence of attractor solution in quintessence model. In this 
paper we take potential Ty(X) as the Mexican hat potential which satisfies the above condition W{Xc) > 0, 
W'(Xc) — and W"{Xc) > 0. Therefore, our quintessence model can admit a late time attractor solution. 

Based on the analysis of the above autonomous system, we can conclude that: Because ^^^^^-^ — jg^ — > 

0, the coupling between matter and dilaton can't change existence of attractor solution in Eqs.(22)(23), that 
is to say, the coupled term can't change the fate of the evolution of the universe. In one word, the existence 

— iacTQX — 3iV 

of couple term ^^^^^ — j^^ — between matter and dilaton affects the evolutive process of the universe, but 

not the fate of the universe. 

In what follows, wc shall show mathematically the evolutions of the components of cosmic density il, the 
evolution of the parameter of state equation cu and the evolution oi X, Y with respect to N in quintessence 
model. The attractor behavior is also shown in phase plane. 

3. Phantom Model 

When we regard dilaton as phantom field, its kinetic energy will be negative. So, the phantom field equations 
become 

H', = l[-l'y' + Wia)+e-'^y] (31) 

dW _ 1 
'da ^ 2 
Based on the same analysis like the quintessence model, we have 



ff + 3Hp<7--j- = -ae-°"'{p-3p) (32) 



^2 ^ h2[ h^' + Wia) ^ ^^^^^_ i„.(^^3 ^ f^^^^(- )4] (33) 

Pc,i a a 



a + 3H& - -— = -ariie-^°"'a-^ 
da 2 

According to the transition (20), the phantom equation fields become 

dX aoY 



dY 



-3Y 



dN H 



dN (7'^H AH 

where the prime denotes the derivative with respect to X, and H can be expressed as 



H = H, 



-\atY^ + W{X) ^ ^^^^^-i..„x^-3A. ^ 0,,a-^^]^ 



The equation of state parameter becomes 



at,Y^ 



-W{X) 



<T*y2 



W{X) 



(34) 

(35) 
(36) 

(37) 
(38) 



Eqs.(35)(36) and their initial conditions determine the evolution of phantom universe and the behaviors 
of the late time de Sitter attractor, which include the existence and the stability of the solution. According 
to the similar analysis to the quintessence model, we also conclude: W{Xc) > 0, W'{Xc) = and W"{Xc) < 
0(a;i < 0, a;2 < 0) is the sufficient condition for existence of attractor solution in phantom model. In this 
paper we take potential W^(A^) as the Mexican hat potential which satisfies the above condition W{Xc) > 0, 
W'{Xc) = and W"{Xc) < 0. Therefore, our phantom model can admit a late time attractor solution. 

In what follows, we shall show mathematically the evolutions of the components of cosmic density fl, 
the evolution of the parameter of state equation lo and the evolution of X, y with respect to N in phantom 
model. The attractor behavior is also shown in phase plane. 



4. Mexican Hat Potential ^(ct^ - e 



Wn 



For this type of Mexican hat potential, it has two extremum points in the range cr > 0: a minimum at 
a — e and a maximum at ct = 0. The non-conventional parameter Wq in this potential, moves the potential 
up and down, which is equivalent to adding a cosmological constant to the usual Mexican hat potential. We 
show the feature of Mexican hat potential mathematically in Fig.l. 



W(cS 




Fig.l Mexican hat potential when we set Wo = 0{lower line), 0.1{middle line), 0.2(upper line). 

The critical point of the above autonomous system is {Xc, 0), where Xc is defined by W'{Xc) — 0. Note 
that the energy density of the dilaton field at the critical point is W(Xc) and should not vanish, thus the 
sufficient condition for the existence of a viable cosmological model with a late time de Sitter attractor 
solution should be that: the potential of the field has non-vanishing minimum value in quintessence model, 
on the contrary, the potential should have non-vanishing maximum value in phantom model. The extremum 
difference of potential between quintessence model and phantom model results from the unusual physical 
feature of phantom field. 

Since the Mexican hat potential(Wo =/= 0) has a non-vanishing minimum value and a non-vanishing 
maximum value, there must exist late time attractors in the quintessence model and phantom model, which 
both drive from dilaton field. Based on this, we study the attractor behaviors of quintessence and phantom 
model in Mexican hat potential. 

In quintessence model, the critical point is Xd = e, and at this point the potential has a minimum value 
Wq. Obviously, it is an attractor solution, which corresponds to an equation of state parameter to = —1 and 
a cosmic density parameter Q„ = 1. The evolution of the parameter of state equation ui, the cosmic density 



parameter fi, Y with respect to N are shown numerically in Figs. 2-4. We also show the attractor property 
for quintessence in phase plane in Fig. 5. 



n. 



N 

Fig. 2 The evolution of cosmic density param- 
eters Q,c{solid line), ilrnidash — dot line), 
Qr{dot line) with respect to N in Mexican hat 



potential for quintessence. We set 
a = 0.01, Wo = 1.00. 
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10" 




^ Fig.4 

The evolution of Y with respect to N for dif- 
ferent coordinate initial conditions Xi = 0.5, 
Yi = 0.5; Xi = 1.2, Yi = 1.2 in Mexican hat 
potential for quintessence. 




Fig.3 

The evolution of the parameter of state equa- 
tion uj with respect to N in Mexican hat po- 



tential for quintessence. 
a = 0.01, Wo = 1.00. 



We set 



Hi 



= 10" 




Fig. 5 The phase portrait for quintessence in 
Mexican hat potential for several different ini- 
tial conditions. The critical point Xc ~ 3.333 
corresponds to a de Sitter solution which is a 
stable spiral 



In phantom model, according to the sufficient condition of existence of attractor solution, we can obtain 

4 

the critical point is Xc2 = 0, and at this point the potential has a maximum value ^^ -I- Wq- Obviously, it is 
also an attractor solution, which corresponds to an equation of state parameter tu = —1 and a cosmic density 
parameter fla- = 1. The evolution of the parameter of state equation lu and the cosmic density parameter 
fi with respect to A'' are shown numerically in Fig. 6 and Fig. 7. We also show the attractor property for 
phantom in phase plane in Fig. 8. 




Fig. 6 The evolution of cosmic density param- 
eters Qcr {solid line), flm{dash — dot line), 
Qr {dot line) with respect to N in Mexican hat 
potential for phantom. We set -S- ~ 10~ 
a = 0.01, Wo = 1.00. 
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N 

Fig. 7 The evolution of the parameter of state 
equation uj with respect to N in Mexican hat 
potential for phantom. We set -^ = 10~^^, 
a = 0.01, Wo = 1.00. 
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Fig. 8 The phase portrait for phantom in Mexican hat potential for several different initial conditions. The critical 
point Xc = corresponds to a de Sitter solution which is a stable spiral 



In what follows, we shall show the effect of nonminimum coupling term 



ar/iC 



AH 



on the evolution 



of X, Y and the attractor solution. From Figs. 9-11 and the analysis in Section 2, we obtain that the existence 
of couple term between matter and dilaton affects the evolution course of the universe but not the result. 
That is to say, either quintessence or phantom model admits an attractor in spite of the existence of coupling 
term. 
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Fig. 9 The evolutive comparison of X with 
respect to A^ in quintessence with coupling 
term(real line) and without coupling term(dot 
line) . 




Fig. 10 The evolutive comparison of Y with 
respect to A^ in quintessence with coupling 
term(real line) and without coupling term(dot 
line) . 
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Fig. 11 The phase portrait comparison for quintessence with coupling term(real line) and without coupling term. 

Both of them admit the same one attractor Xc — 3.33. 



5. Conclusions 

In this paper, we consider dilaton field as quintessence and phantom in dilatonic dark energy model. 
We investigate the existence of a late time attractor in the the two cases. The sufficient condition of the 
existence of attractor solution is: the potential of the field has non- vanishing minimum value in quintessence 
model, on the contrary, the potential should have non-vanishing maximum value in phantom model. The 
physical reason of the condition difference between quintessence and phantom is clarified as follows: For a 
usual dynamic system(quintessence), the system will be stable when the total energy of the system arrives 
at the minimum value. So, for quintessence the sufficient condition is that the potential of the field has 
non- vanishing minimum value. However, for a phantom system whose kinetic energy is negative, the system 
will be stable when the total energy of the system arrives at the maximum value. So, for phantom the 
sufficient condition is that the potential of the field has non-vanishing maximum value. 

The Mexican hat potential has a local minimum and a local maximum which makes us consider quintessence 
and phantom in the same potential. In quintessence model, the critical point is Xd — e, and at this point 



the potential has a minimum value Wo- In phantom model, we obtain that the critical point is Xc2 — 0, and 

4 

at this point the potential has a maximum value ^^ + Wq. We show mathematically that both quintessence 
model and phantom model admits a late time attractor solution corresponding to an equation of state 
Lj = — l(Fig.3 and Fig. 7) and a cosmic density parameter H.^ — l(Fig.2 and Fig. 6), which are important 
features for a dark energy model that can meet the current observations. Such evolution behaviors will avoid 
the cosmic doomsday. Since Teq ~ 5M{Qoh'^)eV ~ 2.843 x lO^K, To ~ 2.7K, ai = 1, the scale factor at 
the present epoch ao would nearly be 1.053 x 10'', then we know iVo = Inao — 9.262. According to No, we 
obtain the current value ^a,quintessence — 0.707036 and fla-,phantom — 0.703174, which both meet the current 
observations well. 

-^a^QX -3JV 

Finally, we analyze the effect of the coupling term ^^^^^^ — jjj — on the evolution oi X , Y and 

attractor(Figs.9-ll). We find that it changes the evolutive process of X and Y, but not the existence 
of the attractor solution. These mathematical results show that, because of the existence of this coupling 
between dilaton and matter, the evolution of X and Y with respect to N in quintessence model will be faster 
than those of the other models without coupling. According to the same analysis methods, we can expect 
that the coupling term affects the evolution of w and Q in the same way. 

We can find that the state parameter of dilaton as quintessence and phantom can both evolve to -1 from 
values greater than -1 and from values smaller than -1, respectively. Can it cross over -1? Or dilaton can be 
a candidate of quintom[31]? Vikman has proved that the crossing of w = —1 is classically forbidden in the 
single scalar field models[32]. However, in the work of Singh [33], we can observe such a behavior of crossing 
of ui = — 1 in semi-classical Loop Quantum Cosmology very naturally. If one considers the quantum effects, 
such a crossing is possible for a single scalar field. 
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